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ABSTRACT 



Context. Symbiotic X-ray binaries (SyXBs) are a recently discovered subclass of low mass X-ray binaries. Their growing number 
makes them an important evolutionary channel of X-ray binaries. 

Aims. Our goal is to perform spectral analysis and classification of the proposed counterparts to IGR J16358-4726 and IGR J16393- 
4643 and establish their nature as X-ray systems. 

Methods. We used the ESO/UT1 ISAAC spectrograph to observe the proposed counterparts to the two sources, obtaining A"-band 
medium resolution spectra (R = 500) with a S/N > 140. Data reduction was performed with the standard procedure. We classified 
them by means of comparison with published atlases. We performed SED fitting to refine the spectral classification. 
Results. The two counterparts clearly exhibit the typical features of late-type stars, notably strong CO absorption bands in the red part 
of the spectrum. With information from previous X-ray studies, we classify the two systems as two new members of the SyXB class. 
For IGR J 16393-4643, we considered the most probable counterpart to the system, although three other objects cannot be completely 
discarded. For this system, we compared our findings with available orbital solutions, constraining the orbital parameters and the mass 
of the companion star. 

Conclusions. By including two more systems, we increased to eight the number of known SyXBs, which emerges as a non-negligible 
category of galactic X-ray binaries. 

Key words. X-rays: binaries — Stars: neutron — Stars: individual: IGR J16358-4726; IGR J16393-4643 



1. Introduction 

Low-mass X-ray binaries (LMXBs) are interacting systems con- 
sisting of an accreting compact object and a low-mass (typically 
< 1 Mq) main-sequence or slightly evolved late-type star. A 
new sub-class of LMXBs has emerged, in which the secondary 
is an M-type giant. These ra re systems are calle d symbiotic 
X-ray binaries (SyXBs) after Maset ti et alj d2006l) by analogy 
with symbiotic binaries, in which a white dwarf accretes matter 
from an M giant companion, either from the wind of the M star 
or via Roche-lobe overflow. At present, only si x X-ray sources 
are known to belong to this subcla ss: GX 1+4 (Davidse n et all 
1 1977 1 : IChakrabartv & Roche] 1 19971) . 4U 1700+24 (iGarcfa et alj 
1983[ IMasetti et all 120021) . 4U 1954 +319 dMasetti et all 
200(1 lMattanaetalJ l2006l)7 Scutum X -l dKaplan et all 120071) . 
IGR J16 194-2810 dMasetti etal1l2007l). and 1 RXS J1 8043 1.1- 
273932 dNucita et all 120071) . IMasetti et alj d2007l Table 3) 
present ed a summa r y of the main properties of the first 5 
SvXBs. lNucita et all d2007l) including the last discovered one. 

The increasing number of systems with an evolved giant 
donor makes this subclass an emerging evolutionary channel of 
X-ray binaries. They have also been prop osed as the probable 
proge nitors of most wide-orbit LMXBs (Chakrabart v & Rochd 
fl997h . 



The six-years of INTEGRAL/IBIS data are improving 
substantially our knowledge about G alactic X-ray binaries (see, 
for instance, iBird et a II 120071 120101) bv rdetecting a new class 
of heavily absorbed supergiant m assive X-ray binaries (first 
proposed bv iRevnivtsev et al1l2003l) . by detecting and allowing 
the stud y of supergiant fast X-ray transients (SFXTs; see, for 
instance, Smith et al. 2006; Sguera et al. 2006; Negueruela et al. 
120071 ISidoli et al.1 120071 iNesnoli et alJ l2008d) . bv doubling 
the number of kno wn high-mass X-ray binaries (HMXBs, see 
Walt er et all 12006), and by discovering a vast number of new 
magnetic cataclysm ic variables (CVs; see iBarlow et all 12006; 
IMasetti et all 12009). Among INTEGRAL r esults, one can als o 
identify one new SyXB, IGR J 16 194-28 10 (Mas etti et alj2007b . 



Send offprint requests to: Elisa Nespoli 

* Based on observations collected at the European Southern 
Observatory, Chile (Programme ID 079.D-0668-B) 



IG R J16358-4726 was discovered by IRevnivtsev et all 
(2003), then serendipitously observed by Chandra, which lo- 
cated the source at RA = 16 *35 m 53.8 J , Dec = -47°25 '41.1" 
with an accuracy of 0.6" dKouveliotou et al] 120031) . This 
allowed these authors to identify an infrared counterpart, 
2 MASS J16355369-47 25398. Subsequent astrometry performed 
by IChatv et all d2008l) . although stating that the 2MASS ob- 
ject is 1.2" from the Chandra detection, resulted in establish- 
ing it as the true counterpart, mainly because of its brightness 
(J = 15.41 mag, H = 13.44 mag, K = 12.59 mag). The X-ray 
source was identified as a transient, and its spectrum fitted with 
a power law of spectral index T ~ 0.5, hydrogen column den- 
sity Nu_ 1 33 x 10 23 cm , and a FeKa fluorescence line at 6.4 
keV dPatel et alJl2007l) . From accurate spectral and timing anal- 
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Table 1. ISAAC/VLT journal of observations. In the fourth column, we report the net accumulated exposure time. Column five 
gives the obtained signal-to-noise ratio. The references list in the last column relates to the identification of the optical/infrared 
counterpart. 



Source 


K mag 


Start time (UT) 


Exp. time (s) 


S/N 


IR Counterpart 


Reference 


IGR J16358-4726 
IGR 116393-4643 


12.6 
12.8 


2007-04-05 08:26 
2007-04-05 08:56 


480 
1280 


120 
160 


2MASS J16355369-4725398 
2MASS J16390535-4642137 


Kouveliotou et al. (2003) 
Bodaghee et al. (2006) 



ysis of multi-satellite data. lPatel et al.l (120071) also proposed that 
the source is a neutron star, estimating its magnetic field to be 
between 10 13 and 10 15 G, if the retrieved spin up was caused 
by disc accretion, and thus proposing that the object might be a 
magne tar. A spin period of 5880 + 50 s w as found by Patel et all 
( 200 4]) and confirmed b y several authors dLutovinov et al.| [2005; 
iMererfietti et~ai1 12006)." The source is understood to be located 
in the Norma Galactic arm, at a distance of 5-6 kpc or 12-13 
kpc, depending on which crossing of the arm with the line of 

sight is chosen dLutovinov et al.ll2005h. 

From NIR spectroscopv JChatv et al.l (12008) inferred that the 
system is a HMXB, and proposed that the counterpart might be 
a sgB[e] star because of the detection of H Brackett series, He I, 
and Hell absorption lines and a forbidde n [Fell] line at 2.22 
pm. This result was then corroborated by Rah oui et al.l ((2008) 
by means of SED fitting. 



ISugizaki et all J200lb first discovered IGR J16393-4643 
with the ASCA sateUite and listed it as AX J16390.4-4642. 
During the first I NTEGRAL scan of the Galactic plane, it 
was later detected (M alizia et al.l 12004b . Subsequent observa- 
tions with XMM-Newton/EPIC derived an improved position 
(RA = 16 A 39 m 05.4 J DEC = -46°42T2"), which enabled 
iBodaehee et all J20061) to identify 2MASS J16390535-4642137 
as the most probable counterpart to the source. These authors fit- 
ted the X-ray spectrum with a highly absorbed (Na = 2.5 x 10 23 
cm -2 ) power law (F = 0.8 + 0.2), which, with the discovery of a 
pulse period of 912.0 + 0.1 s, led to co nclude that the source is 
an X-ray pulsar dBo daghee et al. 2006). 

iThompson eTal] d2006l hereafter T06) used RXTE data 
to carry out a pulse timing analysis to determine the orbital 
parameters of the system. They obtained three mathematical 
solutions, among which the most plausible one indicated, with 
an orbital period of 3.7 days and a mass function of 6.5 ±1 Mq, 
that the system should be a HMXB. The other two solutions 
foresaw the following parameters: an orbital period of 50.2 
days for a mass function of 0.092 Mq and an orbital period 
of 8.1 days for mass function of 221 Mq; both of them were 
rejected by the aut h ors as being statistically and physically 
weaker. Chat v et all d2008l) obtained accurate astrometry and 
photometry of the field of IGR J16393-4643, finding three more 
candidate counterparts, although their lower NIR brightness 
seemed to favor the 2MASS object. They also performed SED 
fitting of optical to MIR observations, inferring a BIV-V spectral 
type for the companion. 

In this paper, we present Zf-band spectroscopy of the pro- 
posed counterparts to the two X-ray systems. In the case of 
IGR J16393-4643, our data represents the first available IR spec- 
trum; in the case of IGR J16358-4726, our data is the first IR 
spectrum beyond 2.3 /mi. We demonstrate the need for a new 
classification of the counterparts, which implies a new classifi- 
cation of the systems as SyXBs, if the proposed counterparts are 
the true ones. In the next section, we describe our observations 



and data reduction; in Sect. 3, we report the obtained spectra, 
analyze their features, and propose a classification; in Sect. 4, 
we discuss our results, and in Sect. 5 we point out the conclu- 
sions. Preliminary re sults of our data analysis were published in 
iNespoli etalJd2008al) . 



2. Observations and data analysis 

The two proposed counterparts to the X-ray sources were 
observed in s ervice mode on 2007 April 5 with the ISAAC 
spectrograph dMoorwood et all 1 19981) on UT1 at ESO/Paranal 
observatory. The sky was clear during the observations, the 
seeing was < 1 .4", and the targets were observed at airmasses of 
1.08 and 1.09, respectively. Data were taken in the short wave- 
length - low resolution mode with a pixel scale of 0.147"/pixel 
and a resolution of 500. Table [TJ reports the observation log, 
including the achieved signal-to-noise ratio (S/N). Typical 
on-source integration times for standard stars were between 6 
and 10 seconds. 

Data reduction was performed using the IRAB3 package, fol- 
lowing standard procedures for IR spectra. We first corrected 
for the inter-quadrant row cross-talk, a feature that affects the 
ISAAC detector. We then performed sky subtraction, applied 
dome flat-fields, extracted and rectified the one dimensional 
spectra. Wavelength calibration was accomplished using xenon 
and argon lamp spectra. Spurious features, such as cosmic rays 
or bad pixels, were removed by interpolation, when necessary. 
The reduced spectra were normalized by dividing them by a fit- 
ted polynomial continuum. 

To ensure accurate removal of atmospheric features from 
th e spectra, we followed a strategy similar to that descripted 
bv lClark & Steelel {2000). Both an AO - A3 III-V and a G2-3 V 
standard stars were observed immediately before or after each 
target to obtain very small differences in airmass (differ- 
ences of 0.008 and 0.007 in airmasses were accomplished for 
IGR J16358-4726 and IGR J16393-4643, respectively). To com- 
pute the telluric features in the region of the HI 21661 A 
(Brackett-y line, or Bry), which is the only non-telluric feature 
in the A-star spectra, we employed the observed G-star spectra 
divided by a solar spectrumj that had been properly degraded in 
resolution and dispersion corrected with the dispersion solution 
obtained for the ISAAC spectra. The spectra of the A star, G star, 
and the solar one were then aligned in wavelength space. A tel- 
luric spectrum for each science target was obtained by patching 
into the A-star spectrum the ratio of the G star to the solar spec- 
trum in the Bry region (we selected the range 21 590 - 21 739 



1 IRAF is distributed by the National Optical Astronomy 
Observatories which is operated by the Association of Universities for 
Research in Astronomy, Inc. under contract with the National Science 
Foundation 

2 We used the NSO/Kitt Peak FTS solar spectrum, produced by 
NSF/NOAO 
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Fig.l. £-band spectrum of 2MASS J16355369-4725398 and 2MASS J16390535-4642137, the infrared counterparts to 
IGR J16358-47262 and IGR J16393-4643, respectively (first two spectra from above). The positions of identified spectral fea- 
tures are marked by solid lines. With the aim of comparison, we also show the /f-band spectrum of a blue supergiant (third spectrum 
from above). The three targets were observed in the same night with the same telescope/instrument configuration. Data reduction 
was performed in the same way. 12 CO absorption bands (A > 22 900 A), typical of late-type stars, are completely absent in the third 
spectrum. 



A). We then corrected for telluric absorption by dividing each 
science spectrum by its corresponding telluric spectrum. A scale 
and a shift factor were applied to the calibration telluric spec- 
trum, to correct in the most accurate way for the airmass dif- 
ference and the possible constant wavelength shift; the optimal 
values for these parameters were then obtained using an iterative 
procedure that minimizes the residual noise. 

3. Results 

We present the spectral analysis and classification results for 
the two targets. The field of NIR spectral classification has not 
yet been fully developed and the level of S/N and resolution 
required to perform a quantitative profile analysis are very 
high (R ~ 12 000 and S/N > 250). Our analysis is therefore 
qualitative, based on comparison with available NIR spectral 
atlases. 

The spectra are shown in Fig. Q] (first two spectra from 
above), while their most remarkable features are reported in 
Table [2] Errors in the equivalent widths are evaluated as the 
deviations from their mean values, measured using different 
values for the continuum. For the sake of comparison, we 
included in Fig. Q] the spectrum of IGR J16493-4348 (third 
spectrum from above), a supergiant X-ray binary (SGXRB) 
that we observed during the same night and conditions as the 
previous ones. The two first spectra are very similar, their most 
prominent features being the series of strong CO absorption 
bands between 2.29 and 2.40 fan, which are characteristic of 
late-type stars. Several metallic lines are clearly recognizable, 
although at our resolution they are generally blended. Although 



weaker than the CO bands, they appear significant within the 
measured errors. 

We also tentatively detected He I 20581 A and Bry in faint 
emission in both spectra; since they are not a typical spectral 
feature of late-type stars, these lines might be the signature of 
an accretion disk around the compact object. However, they are 
very weak, and we cannot exclude that they might be a residual 
artifact from the telluric lines removal procedure, which for this 
region should be performed with great care. This interpretation 
is consistent with the almost identical profile of the Bry line 
in the two spectra. In any case, emission features typical of 
X-ray binaries do not seem to be present in almost all of the 
optical spectra of SyXBs (see iMasetti et alj|2006l 120071) . This 
is because of the high total luminosity of these evolved giant 
companions, that should overwhelm the light derived from the 
reprocessing of X-ray irradiation. 

Our ISAAC spectra clearly infer that both objects are of late- 
type and exclude that the systems are HMXBs with OB com- 
panions. From the identified spectral lines, our conclusion is 
that both stars are of K or M spectral type (see, for instance, 
Kleinmann & Halll 11986b iForster Schreiberi l2000t flvanov et alJ 
120041 for a comparison with atlases of cold stars). The forbid- 
de n [Fe II] transition at 22200 A detected in IGR J16358-4726 
bv IChatv et al.1 (120081) is not present in our spectrum. The gen- 
eral morphology of the spectra and the comparison between 
the relative strength of CO 22 900 A and both Nal and Cal 
lines, which are luminosity indicators (Ra mirez et alJ[l997D . in- 
dicate that both stars are giants. In particular, from the empiri- 
cal relation EW(CO) versus. r e ff that these authors found, our 
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Table 2. K-band line identifications for IGR J16358-4724, de- 
noted by (I), and IGR J16393-4643 (II). 



Feature Transition Wavel. (A) EWi (A) EW„ (A) 



He I 


2s 1 S -2p l P° 


20581 


-1.8+0.4 


- 


Mgl 


ir/..: ?>..'(;': ; , 


21066 


0.8±0.2 


1.6±0.3 


All 


4p 2 P° /2 -5/5 1/2 


21099 


_ 


1.8+0.3 


Bry 


4 2 F" - T-G 


21661 


-2.8±0.2 


-0.9±0.2 


Nal 


4rS 1/2 - 4p 2 P° l/2 


22 100 


1.9+0.4 


2.0±0.8 


Fel 


x 5 F° - e 5 D 2 


22263 


1.6+0.2 


0.8+0.2 


Fel 


x 5 F° - e 5 Z), 


22 387 


1.2±0.4 


1.1 ±0.4 


Cal 


4d 3 D 3X1 -AfFl 3 2 


22 636 


1.4+0.4 


1.1 ±0.4 


Mgl 


4d 3 D 3 , 2 , 1 -6/ 3 F° 34 


22 800 


2.7+0.5 


1.6±0.4 


12 C0 


(2, 0) bandhead 


22 900 


9.9±0.3 


11.9±0.4 


12 co 


(3, 1) bandhead 


23 226 


13.0+0.6 


12.9±0.6 


13 co 


(2, 0) bandhead 


23 448 


2.5±0.4 


2.5±0.8 


12 co 


(4, 2) bandhead 


23 524 


7.8±0.4 


9.8±0.3 


13 co 


(3, 1) bandhead 


23 739 


4.2+0.2 


4.5±0.8 


12 co 


(5, 3) bandhead 


23 829 


12.8±0.6 


8.05±0.7 


13 co 


(4, 2) bandhead 


24037 


3.6±0.5 


1.98±0.3 


12 co 


(6, 4) bandhead 


24156 


10.1+0.8 


7.0±0.4 



EW( 12 CO(2, 0)) infers a K giant classification for both objects, 
and nominal values of r eff ~ 4400 K for IGR J16358-4726 and 
T eff ~ 4250 K for IGR J16393-4643. 

We believe that the level of precision achievable in our 
spectral classification cannot be higher than this, because of the 
lack of reliable T e g indicators in the K-bwA in this spectral type 
range. Adequat e indicators in the N IR domain are present in the 
J and H bands (Wright et al. 2008). Unfortunately, we have no 
spectra for these ranges. 

To corroborate the accuracy of our analysis and to avoid 
any possible mistaking of CO absorption bands with telluric 
components, we included in Fig. [TJ besides the spectra from 
the two sources described in this work, the spectrum of 
IGR J 16493-4348, a SGXRB that we observed in the same 
night and conditions as the previous ones. We classified the 
count erpart to IGR J16493-4348 as a B0.5Ib star dNespoli et alJ 
2008bj) and, as one can see, the telluric corrected spectrum of 
the source shows no evidence of the deep absorption features 
that characterize the spectra that are the subject of this work. 

In Figs. [2] and [3] we present the spectral energy distribution 
of IGR J16358-4724 and IGR J16393-4643, respectively, in the 
lo gf/iFj (A) I AkFi (K)) - lo g A plane, us i ng ph otometry given 
by IChatv et all (12008b and iRahoui etail d2008l) . This includes 
photometry in the bands RIJHK and the GLIMPSE 3.6 fim 
and 4.5 jt/m bands for both sources; for IGR J16393-4643, V 
photometry was also available. We built SEDs corresponding to 
different spect ral types, according to the intrin sic colors g i ven by 
iPicklesi (119981) for optical/NIR bands and bv lKoornneel (1 1983k 
for MIR. In the figures, data are marked by crosses, and fitted 
SEDs by lines. Photometric data and spectral type distributions 
are normalized in the K band. We have did not attempt formal 
fitting procedures. The aim of the figures is to illustrate that 
very different choices of spectral type and reddening produce 
reasonable fits to the data. In o ur calculation, we emplo yed the 
mean extinction law, Ry = 3.1 dRieke & Lebofskvfr9 85). 

We plot a blue supergiant and two red giants as represen- 
tative examples of possible spectral types. We found that, even 
with slight differences, every spectral type can be fitted to the 
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Fig. 2. Spectral energy distribution of IGR J16358-4724 
(crosses) and SED obtained from intrinsic colors for three dif- 
ferent spectral types (lines). The errors for the photometric data 
are within the size of the marks. 

data by choosing the appropriate value of extinction. In the 
case of IGR J16358-4726, we obtained best fits by choosing 
A v - 20 mag for a B0I type, A v = 14 mag for a K0III type, 
and Ay =13 mag for a M0III type. Analogous results were ob- 
tained for IGR J 16393-4643 with A v = 12 for a B0I, A v = 8 for 
a K0III and A v = 6 for a M0III. 
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Fig.3. The same as Fig. El for IGR J16393-4643. The errors in 
the photometric data are within the size of the plot symbols. 

The reason why very different spectral types are able to 
reasonably fit the data is the degeneracy between the direction 
of the reddening vector and the effective temperature variation 
in the spectral range that we employed, which makes almost 
every spectral type possible if the opportune extinction value is 
chosen. The results provided by SED fitting do not therefore 
allow the spectral classification to be improved, but provide a 
valuable tool with which to calculate the interstellar extinction 
of any assumed spectral type. 

In the case of IGR J16393-4643, the low values obtained 
for the extinction are inconsistent with the high hydro- 
gen column density measured in the X-ray range, namely 
Nu = 2.5 x 10 23 cirT 2 . This may indicate that the star is not 
the true counterpart to the system; another explanation could 
be an additional source of extinction intrinsic to the system, 
which only affects the compact object where the X-ray emission 
originates. Enormous amounts of absorption, with A^h ~ 10 23 
cm" 2 , has been observed for other SyXBs, (see, for instance, 
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iMattana et afl l2006). 

For the aforementioned spectral types and corresponding 
extinction values, we calculated the distance using the relation 
Mk — K + 5 - 5 log d - Afc. In each case, absolute magn itude s 
and intrinsic color indices were obtained from Pickles 0998). 
The estimated distance values were, for IGR J16358-4726: 36.6 
kpc (if BOI), 1.4 kpc (if KOIII), and 20.1 kpc (if MOIII); for 
IGRJ16393-4643:24.1kpc (if BOI), 1.0 kpc (if KOIII), and 13.9 
kpc (if MOIII). For both sources, a type BOI counterpart would 
thus provide an unrealistically large estimate of the distance; the 
same occurs, for IGR J16358-4726 only, in the case of a MOIII 
companion. 

The remaining estimates are compatible with both the 
classical (four arms, Vallee 2008 ) and the mo re recently pro- 
posed (two main arms, Churchw ell et al.l 1 20 09) spiral structure 
models of our Galaxy. A KOI companion would locate both 
sources in the Sagittarius/Carina inner arm according to the 
classical model, or the Sagittarius minor arm, according to the 
most recently proposed model. Finally, a MOIII companion for 
IGR J16393-4643 would imply that the source is placed in the 
tangential region of the Norma arm with both models. 

The error estimation for the distance depends mainly on 
the error in the mean absolute magnitude assigned to the 
spectral type, which can be very large and reach more than 50% 
when propagated to the distance calculation. Since the derived 
extinction value for each type would also apply to a wide range 
of subtypes, the consequent distance value range would largely 
extend. The distance inferred should thus be interpreted with 
care. 



4. Discussion 

On the basis of the spectral classification of the proposed 
counterparts as cold giants/supergiants and the X-ray features 
of the systems illustrated in Sect. 1, we propose that both of 
them belong to the class of SyXBs, in which the companion to 
the neutron star is an evolved late-type star. In the discussion 
that follows, we use the terms giant and supergiant in its 
physical sense, i.e., related to the mass, internal structure, and 
evolutionary status of the star and not to its luminosity class. 
As discussed above, a precise determination of the spectral type 
and luminosity class is difficult using our K-band spectra. 

A red giant represents the late evolutionary stage of either a 
low-mass star (0.8 to 2.3 Mq) or an intermediate-mass star (2.3 
to 10 Mq). The red giant phase includes the shell hydrogen- 
burning phase (first giant branch), the core helium-burning 
phase (horizontal branch or "red clump" for low-mass stars, 
"blue loop" for intermediate-mass stars), and the shell helium- 
burning phase (asy mptotic giant branch). W e refer, for instance, 
to llbenl (119741) and llben & Renzinil (1 19831) . Red supergiants are 
short-lived stages in t he late evolution of moderately ma ssive 
stars (10 to 25 Mq, see lMassev et alj|2009l: TFiger et al.ll2006l) 

In the case of IGR J16358-4724, we classified the proposed 
counterpart as a K-M giant/supergiant. Strong support for this 
system hosting a neutron star is provided by the detection of 
X-ray pulsations dPatel et al.1 12004) in addition to the description 
of the X-ray energy spectrum as an absorbed power la w model, a 
high energy cut-off, and an Fe line (Patel et al. 2007), which are 



typical features of neutron stars in Galactic X-ray binaries. The 
designation of 2MASS J16355369-4 725398 as a counterpart 
was confirmed in the literature (see iKouveliotou et alj [2003: 
IChatv et al. 2008). Previous NIR spectroscopy from Chatv et al. 
(2008) was used to classify it as a sgB[e], but, as the authors 
warned, their NIR spectrum was very faint and extremely 
absorbed. By ending at 2.3 yum, it also missed the last and 
most significant portion of the K band, where CO bands clearly 
indicate the K-M spectral type of the star. We propose that these 
factors could lead to a misclassification of the system. 

In the case of IGR J16393-4643, our spectrum is also 
indicative of a K-M giant/supergiant. The presence of a pulsar 
in the binary system is validated by both th e spectral shape of 
the X-ray emission dBodaghee et alj [2006) 

and the detection 

of pulsations (T06). Because of its magnitude, the optical 
counterpart that we considered is the most probable one among 
the four objects within the error circle of the X-ray detection. 

According to T06, IGR J 16393-4643 has an orbital period 
of 3.7 days and a mass function of 6.5 Mq (orbital solution 1 
in T06). It seems unlikely that an intermediate-mass red giant 
can fit in such a narrow orbit. To check this, we use the physical 
parameters of a 7 M q>, star along its evolutionary track given by 
ISchalleret all (fl992h . 

The size of the Roche lobe of a star in a binary system is 
usually characterized by the Roche radius R^, defined as the ra- 
dius of a sphere that has the same volu me as the Roche lobe. To 
calculate the Roche radius, we used fhe lEggletonld 19831) approx- 
imation 



R, 



0.49 q 



2/3 



0.6 q 2 ' 3 + ln(l + q 1 ' 3 ) 



(D 



where a is the semimajor axis and q is the mass ratio of the 
binary components (q = M c /Mx). We used the common approx- 
imation for the mass of a neutron star, Mx = 1.4 Mq. 

To obtain the Roche radius for solution 1 of T06, we first 
calculated the major axis of the orbit from the third Kepler law 



An 2 a 3 



M x + M c = — — - , 

G p l 



(2) 



where G is the gravitational constant that we wrote as 
G — 2945.02 Rq Mq d 2 , a is the semimajor axis, and 
p = 3.7 d is the orbital period. For the above parameters, we ob- 
tained a Roche ra dius of Rl - 10.68 Rq . 

According to ISchaller et"aT1 (I 1992b . the radius of a 7 Mq 
star at the beginning of the giant branch is 86.5 Rq. During 
its later evolution, at the core He-burning phase, the minimum 
radius reached is 50 Rq. The only evolutionary phase in which 
a 7 Mq star can fit wihin a 10.68 Rq Roche lobe is during the 
main sequence, for temperatures higher than 15 000 K. This 
corresponds to a B-type spectrum, which is completely rejected 
by our spectral classification. This result shows that the system 
IGR J16393-4643, in the binary orbit proposed by T06, cannot 
harbor a late-type giant star of ~7 Mq. 

Among the orbital solutions calculated by T06, the above 
discussion thus excludes solution 1. Our conclusion may be 
supported by the non-detection of any orbital period of the 
order of a few days (as solution 1 estimates) by Bodaghe e et al.l 
(2006) from ISGRI light curves. We also exclude solution 3 
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Table 3. Results from the analysis of the orbital parameters of 
IGR J16393-4643 for different values of mass of the companion 
star M c . Tabulated values for and L are taken from Schaller et 
al. 1992, for a star at the beginning of the giant branch. For each 
case, the last column shows the corresponding spectral type. 



Table 4. Synoptic table reporting the main known parameters 
for IGR J16358-4724 and IGR J16393-4643. The X-ray lumi- 
nosity is given in the 2-10 keV band. The mass accretion rate 
M was estimated assuming typical parameters for a neutron star, 
i.e. R x = 10 km and M x = 1.4 Mq. 



(M Q ) 



(*©) 



log T eB log(L/L©) 



Re 

(Re) 



Sp. t. 



Parameter 



IGR J 16358-4724 



IGRJ16393-4643 



1 


76.5 


26.8 


3.72 


0.33 


1.75 


54.0 


G5 III 


2 


85.9 


35.2 


3.71 


1.23 


5.19 


30.7 


G5 III 


3 


93.6 


41.7 


3.70 


1.95 


12.46 


23.8 


G8 III 


4 


100.3 


47.3 


3.68 


2.46 


24.35 


20.3 


K0III 


5 


106.1 


52.1 


3.67 


2.85 


40.93 


18.1 


Kl III 



because of, as the authors point out, its statistical weakness 
and physically unrealistic constraints. Although statistically 
less favored than the first, solution 2 appears to be the only 
acceptable solution. With a mass function of 0.092 Mq and an 
orbital period of 50.2 days, it is indeed compatible with our 
spectral classification, as we now show. 

To constrain the mass of the companion, we calculated, 
for increasing values of mass, the corresponding radius at 
the beginning of the giant branch, and compared it with the 
maximum permitted one, i.e., the Roche radius. 



Spectral type 
of the secondary 

V-band magnitude 

A v (mag) 

Distance (kpc) 

X-ray spectrum 

L x (erg s" 1 ) 

I ^secondary 

M(g s" 1 ) 

P spin y0 

P spin I P spin (S ) 

Po* (d) 



K-MIH[1] 
> 23.67 [2] 
12-13 [1] 
5-6 ; 12-13 [31 
[5] 

3 x 10 32 - 2 x 20 36 [5] 

1.6 x 10 12 — 1.1 x 10 16 [1] 
5850 [6] 
3.1 xlO" 8 [5] 



K-MIII[1] 
21.53 [2] 
4-6 [1] 
~ 10 [4] 
[4] 

2 x 10 34 [4,11 

1.1 x 10 14 [11 
-912 [4,7] 
1.0 xlO"" ? [7,1] 
50.2 ? [7,1] 



References: 

[1] this work; [2] Chaty et al. (2008); [3] Lutovinov et al. (2005); 

[4] Bodaghee at al. (2006); [5] Patel at al. (2007); [6] Patel et al. (2004); 

[7] Thompson et al. (2006). 



We again we first calculated the major axis of the orbit by 
applying the third Kepler law in Eq. (0, using p = 50.2 d from 
solution 2 as the orbital period. 

Taking the common approximation Mx = 1.4 Mq, we 
solved E. [2] with respect to a for different values of M c , and 
employed the retrieved value in Eq. Q] to obtain the corre- 
sponding Roche radius. We then considered a red giant at the 
beginning o f the giant br a nch, t ook tabulated values of T e ff 
and L from ISchaller et alJ dl992h . and solved the black-body 
equation, L = 4nR 2 o-T* s , to obtain the stellar radius. For each 
value of mass, starting with M c = 1 Mq, we compared the 
retrieved radius with the Roche radius of the system. Results 
from this analysis are reported in Table [3] where we also show, 
for each case, the corresponding spectral type, obtained from 
temperature and luminosity, and the inclination i of the system 
as calculated from the mass function of solution 2. 

We emphasize that the parameters R c and i in Table [3] 
correspond to the beginning of the first giant branch ascent. In 
consequence, they only represent a lower limit to the radius and 
upper limit to the inclination angle for stars of the given mass. 

Up to M c = 5 Mq, a red giant is able to fit within the orbit 
foreseen by solution 2 of T06, while for higher values the stellar 
radius would exceed the Roche radius of the system. We thus 
set an upper limit of 5 Mq to the mass of the companion star. 
This mass limit excludes the possibility that the companion is a 
supergiant since the lower limit to the mass of a red supergiant 
is 10 Mq as discussed above. 

An alternative explanation is that 2MASS J16390535- 
4642137 is not the true optical/infrared counterpart to 
IGR J16393-4643. A sp ectroscopic study of the other three stars 
(see IChatv et ail 120081 Fig. 2c) found in the error box of the 
X-ray source is needed to confirm or exclude this case. If one 



of those targets were an OB supergiant, this would then support 
the classification of the system as a HMXB, and would make 
the orbital solution 1 from T06 acceptable. 

For the sake of comparison, we collect in Table [4] the main 
pr operties of the two s ystems studied, thus extending Table 3 
bv lMasetti et al.l (120071) . The mass accretion rate was calculated 
from the X-ray luminosity, as 



where R ~ 10 km is the radius of the neutron star, and M ~ 
IAMq its mass. For IGR J16393- 4643, L x was obtained from 
the X-ray 2-10 keV flux given by iBodaghee et alJ (12006). and 
assuming a distance of 1 .4 kpc, according to our estimate in the 
most probable case of a K0III counterpart. 

Although the class of SyXBs is not yet well characterized we 
note that the known properties of these two systems and those of 
the other six members are homogeneous. We confirm the pecu- 
liarity of the class with respect to the spin, which exhibits pe- 
riodicities spanning a wide ra nge, from hundreds of seconds to 
hours (see Mase tti et alJl2007l and references therein). 

5. Conclusions 

We have used infrared medium-resolution spectroscopy to 
classify the proposed counterparts to IGR J16358-4724 and 
IGR J16393-4643. In this work, we have found that: 

- both objects exhibit the typical features of cold giant or su- 
pergiant stars. Our classification constrains the spectral type 
to be K-M. 

- previous classifications of the systems as HMXBs have been 
dismissed. 



E. Nespoli et al.: IGR J16358-4726 and IGR J16393-4643: two new SyXBs 



7 



- spectral classifications of the counterparts as late type-giant 
stars, and an X-ray behavior characteristic of neutron stars 
known from previous works, has allowed us to include both 
systems in the SyXB class, increasing the number of this 
group from six to eight members. 

- am ong the orbital so l utions for IGR J16393-4643 proposed 
by iThompson et al.l (120061) from a pulse-timing analysis, 
only one, with a mass function of 0.092 Mq and orbital pe- 
riod of 50.2 days, is compatible with our spectral classifica- 
tion. 

- an upper limit to the mass of the companion star in 
IGR J16393-4643 was set to be 5 Mq. 

The results that we have obtained for IGR J16393-4643 are 
valid if the candidate counterpart is confirmed by excluding as a 
possible companion one of the other three sources found in the 
error circle of the XMM-Newton position. Future spectroscopic 
studies of these faint sources will be able to confirm without 
any doubt the classification of the system. Our data clearly high- 
light the importance of obtaining A'-band spectroscopy beyond 
2.3 fim, where CO absorption features are located. With spectra 
being limited shortward of this range, which can be featureless 
in particular at low resolution and S/N, classification can be dif- 
ficult to perform and misleading. 
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